Topologically protected mid-gap states induced by impurity in one-dimensional 
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Based on the discovery of the nontrivial topology of one-dimensional superlattices, we show that 
midgap states will emerge in such systems induced by a single on-site impurity. Besides the trivial 
bound state located at the impurity site, these midgap states are localized at the adjacent sides 
of the impurity, carrying the information of the topology and behaving like the effect of the open 
boundary conditions in the limit of a large attractive or repulsive impurity potential. Using this 
feature, the impurity can be used to detect the topology of the superlattice system and to realize 
the adiabatic pumping between the opposite sides of the impurity in cold-atom experiments or in 
photonic crystals. 

PACS numbers: 05.30.Fk, 03.75.Hh, 73.21. Cd 
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I. INTRODUCTION 

In recent years, due to the very fast development of 
the cold atom techniques 1 , quantum simulation of topo- 
logical insulators in cold atomic systems has been at- 
tracting more and more attentions experimentally and 
theoretically^. However, while nearly all studies fo- 
cus on the two-dimensional or three-dimensional sys- 
tems, one-dimensional (ID) topological insulators are 
less known besides the Su-Schrieffer-Heeder model 3 . Re- 
cently, the discovery of the topology of ID superlattices 
adds a new member— Because of the easy realization 
of the superlattices both in cold atoms^ and in photonic 
quasicrystals^, it paves the way of studying the topology 
of ID superlattices, to date including topological phase 
transitions 7 , fractional topological states 8 -, topologi- 
cal insulator o 9 ' 10 , topological Mott insulators^, and even 
topological superconductors^. 

As we know that the detection of the topological edge 
states depends on the open boundary conditions (OBCs) 
of the system^ 3 -, in cold atoms the confining harmonic 
trap blurs the boundary effects, making the edge states 
not obvious^. To demonstrate the nontrivial topology 
of the ID superlattices, a variety of proposals have been 
given, such as detecting the plateaus in the density pro- 
file subjected to a harmonic trap in cold atoms^, and 
the direct observation of the localized edge states using 
waveguides^. In this paper, we find that a single on-site 
impurity will induce midgap states for the superlattices 
and they are localized next to the impurity site in the 
limit of the large impurity potential, behaving like the 
edge states under OBCs. Using this feature, we can eas- 
ily detect these states in the interior of the system instead 
of at the boundaries, and also realize the transportation 
of midgap states between the opposite sides of the im- 
purity. By the technique of the spatially resolved radio- 
frequency spectroscopy (SRRFS)^, the local density of 
states (LDOS) can be measured directly. So the impurity 
supplies a more direct and realistic scheme to study the 
topology of the ID superlattices. 



The paper is organized as follows: After a brief in- 
troduction, in Sec. |TT] we first introduce the midgap 
states induced by a single on-site impurity in ID spin- 
less fermionic superlattices and show the relationship of 
these states to the Z-type topology possessed by such su- 
perlattices. In Sec. IIII1 the behavior of the midgap states 
can be used to show the topology of the superlattices by 
the technique of the SRRFS in cold atoms to detect the 
LDOS in the harmonic trap. Sec. IIVI generalizes the 
discussion to the spin-1/2 case and the relationship of 
these states to the if 2 -type topology. Alternatively, Sec. 
M gives another scheme to demonstrate the topology of 
superlattices through midgap states using photonic crys- 
tals and realizes the adiabatic pumping between opposite 
sides of the impurity. In Sec. IVI| we have a summary. 



II. Z-TYPE TOPOLOGICAL INSULATORS 
WITH A SINGLE ON-SITE IMPURITY 

Firstly, we concentrate on the ID tight-biding super- 
lattices for spinless fermions with L lattice sites, doped 
with a single on-site impurity. The Hamiltonian is as 
follows: 



H = -t> (cJc l+ i+H.c.) + 



E 



Vihi + Fn 



(1) 



where cj (cj) is the creation (annihilation) operator at 
site-z for spinless fermions, and hi — cjcj. t is the 
hopping amplitude and set to be the unit of energy. 
Vi = Vcos(2iraia + 6) is the on-site modulated poten- 
tial with V being the strength, 1/a the periodicity, and 
5 the arbitrary phase shift. F is the potential of the im- 
purity. Without loss of generality, we choose odd number 
of lattice sites and put the impurity at the center of the 
chain (the center site is labeled as the origin) in the fol- 
lowing discussion. The lattice spacing, a, is taken to be 
1, for convenience. 

Without the modulation, there is only one bound 
state localized at the impurity site with its energy lower 
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(F < 0) (Fig. ma)) or higher (F < 0) (Fig. □») 
than that of the scatting states. When the superlattice 
is turned on, besides the trivial bound states, there are 
also some states emerging within the bulk gaps. These 
midgap states are also localized. However, they are not 
localized at the impurity site but on the ones nearest to 
it (Fig. [5]). Because the behaviors of the cases for F > 
and F < are similar except for the energies of the trivial 
bound states, we mainly focus on the F < case below. 
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FIG. 1: The single particle eigenenergies in ascending order 
with (a) the attractive (F = —5) and (b) the repulsive (F — 
5) impurities located at the center of the ID tight-binding 
chain without modulation under PBCs. The insets are density 
profiles of the bound states labeled by blue stars, which are 
localized at the impurity site with energies (a) lower or (b) 
higher than that of the scattering continuous. Here we choose 
t = l,L= 201. 
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FIG. 2: For the attractive (F = —5) (a-d) and the repul- 
sive (F = 5) (e-h) impurity with the periodic modulation, 
a — 1/3, the eigenenergies with the midgap states labeled 
by colored triangles are shown in ascending order. The corre- 
sponding density profiles (b-d, f-h) of the midgap states in the 
same color demonstrate the ways of localization at one side 
of the impurity or the impurity site. The other parameters 
we choose are t = 1, L = 201, V = 0.5, S = n/A under PBCs. 



edge states, when we shift the phase of the superlattice 
the midgap states do not transfer from one edge to the 
other, but from one side of the impurity to the other, 
which just happens at a very small space helpful to the 
experimental detection. So naturally the single on-site 
impurity is a good candidate for detecting the topologi- 
cal properties of the superlattice systems, which has been 
proposed to detect zero-mode Majorana fermions in p- 
wave topological superconductor o 15 ' 17 . 



We have known that without the impurity, this ID 
superlattice system is topologically non-trivial^, which 
is clearly seen when we adiabatically shift the phase, 
5, under OBCs the edge states will connect the near- 
est bulk bands and thus transfer from one end of the 
chain to the other. A Z-typc topological invariant, 
Chern number, can be defined in the (fc, <5)-space un- 
der periodic boundary conditions (PBCs)*, i.e., C — 

Jo*^ Jo™ dS[dkAg — dsAk] with the Berry connec- 
tion A s = i(4>(k,5)\d s \4>(k,5)) (s = S,k), where k is the 
quasimomentum and <p{k, S) is the Bloch wavefunction of 
the Hamiltonian for each S. 

Now if we turn on the impurity potential, the midgap 
states appear in each bulk gap varying with the phase 
shift S (Figs. EJa) and[3Je)), and as the strength of the 
impurity potential, V, increases, these states will become 
more and more localized at the adjacent sites of the im- 
purity, and connect the nearest bulk bands just like the 
edge states under OBCs (Fig. [3X e )) - Different from the 



III. DETECTION OF MIDGAP STATES IN 
COLD ATOMS 

To realize the superlattice Hamiltonian (jTJ) in cold 
atoms^, we first trap an atomic gas in a primitive ID 
optical lattice potential, V\{x) — V\cos 2 (kix), with the 
wave number k\ and the strength V± , then superimpose a 
secondary optical lattice potential, V(x) = Vcos(k 2 x+S), 
with wave number ki as a weak perterbation if V\ 3> V. 
In the framework of the tight-binding approximation, the 
primitive and the secondary optical lattices are mainly 
supplying the hopping term and the on-site term of the 
Hamiltonian (JTJ) , respectively. Here a — kijk\ can be 
tuned by the wave numbers of the two sets of optical 
lattices and 5 is the overall relative phase shift between 
themi^. 

We know that the observation of the edge states is dif- 
ficult in cold atom experiments due to the harmonic trap 



3 



(a) F=-3 



(e) 



F=-10 



2 
1 

w 



iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniiiiiiiiiiiiiiiiiiiii 




UJ 




FIG. 3: The eigenenergies versus 8 for a single impurity of 
a weaker (F = —3) (a-d) and a stronger (F = —10) (e-h) 
on-site potentials under the periodic modulation, a = 1/3. 
The midgap states are labeled by colored dots, which tend to 
connect the nearest bulk bands and to localize at one side of 
the impurity in the limit of a large impurity potential. The 
corresponding density profiles show some specific point of S's 
to demonstrate the transfer from one side to the other. The 
other parameters we choose are t = 1, L = 201, V = 1 under 
PBCs. 



FIG. 4: LDOS for spinless fermions with a single impurity of 
a large on-site potential (F — —10) within a slowly varying 
harmonic trap (Vh = 0.0005). The midgap states in the bulk 
gaps are transferred from one side of the impurity to the other 
by tuning the phase shift, (a) 8 = tt/3, (b) 8 = 2tt/3, (c) 
5 = it. Here we take a = 1/3, t= 1,L= 201, V = 1. 



IV. Za-TYPE TOPOLOGICAL INSULATORS 
WITH A SINGLE ON-SITE IMPURITY 



for confining Bose-Einstein condensates (BECs), which 
weakens the boundary conditions and reduces the parti- 
cles at both wings of the trap. Here we show that the 
impurity can induce midgap states and has the similar ef- 
fect of the OBCs in the uniform systems. These midgap 
states hold the information of the topology of this system 
and behave like the edge states only in a very small space 
as we mentioned before and overcome the shortcomings 
from the harmonic trap. So placing an impurity into the 
center of the harmonic trap can transfer the detection of 
the topological edge states from boundaries of the system 
to the interior, and the localization of the midgap states 
will be much more notable than those when they are at 
edges. 

We calculate the LDOS in Lehmann representation at 
zero temperature, 

p(i,u) = (G\a\n) \ 2 5(uj-uj n o) + \ (G\ c\ \n) \ 2 8{uj+uj n0 )] 

n 

_ (2) 

where w„o = E n — Eg with \G) and \n) being the ground 
state with energy Eg and the ra-th excited state with 
energy E n , respectively. This quantity can be directly 
measured experimentally by SRRFS 16 . Fig. |4] shows the 
LDOS after adding the harmonic trap term, Hh = Vhi 2 
with Vh being the strength of the trap. We can see that 
as we tune the 5, the midgap states are firstly localized at 
one side of the impurity, then merge into the bulk bands 
becoming extended, and finally transfer to the other side 
of the impurity. 



For spin-1/2 fermions, a spin-dependent secondary 
optical lattice is added with the form V(x) — 
Vcos{2 , Kax ± S), which can be implemented by two 
counter-propagating laser beams with linearly polariza- 
tions forming an angle 2< 5 19 ' 20 . So Hamiltonian (p} can 
be extended to the spin- dependent one just regarding 
ci = (cVf,Cij,) T as a two-component spinor for spin- 
up and spin-down fermionic operators and Vj as Vi = 
diag(Vcos(27rm + 5), Vcos(2itai — 5)). From the same 
procedure like the ID spinless superlattice to simulate 
the quantum Hall effect and demonstrate the Z-type 
topology, this generalized ID spin-dependent superlat- 
tice can be used to simulate the quantum spin Hall ef- 
fect and demonstrate the ^2-type topological insulators 9 . 
Likewise, here the spin-dependent Chern number can be 
also defined in the (k, <5)-space under PBCs as C a — 

h C q dk C d8[d k A Sa - d s A ka ] {a =t or |) with the 
spin-dependent Berry connection A sa and the Bloch wave 
function <p(k, S) for fixed (5. Using these spin-dependent 
Chern numbers, a total Zi topological invariant is con- 
structed as v — (C-|- — Cj,)/2, which demonstrates the 
^2-type topology. Furthermore, this formalism is also 
applicable to spin-orbit systems where the total spin does 
not conserve, the only change of these definitions for the 
topological invariants is to replace the fermions by quasi- 
femions which are the mixture of the spin-up and spin- 
down fermions with pseudo-spins. As long as the bulk 
gap does not close, the topology of the system will not 
change. 

When the impurity is embedded, because of the pres- 
ence of the spin-dependent potentials of the superlattice, 
different from the spinless case, the system may possess 
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simultaneously spin-up and spin-down fermions localized 
at one side or both sides of the impurity. In principle, we 
can always translate the whole superlattice by a phase, 
5q, which makes the midstates shift wholly by So, i.e., 
Vi = diag(V co s(2-jrai — So + S), Vcos(2irai — Sq — S)). By 
tuning the whole phase So and the relative phase S, we 
can manipulate the emergence of the midgap states at 
one side or both sides of the impurity simultaneously or 
non-simultaneously. For example, without the phase dif- 
ference of the two component fermions, the whole spec- 
trum is of double degeneracy with spin-up and spin-down 
fermions having the same energy, and the midgap states 
will always be localized at the same side of the impu- 
rity for fixed Sq. Another case is that we can always 
make the whole spectrum mirror-symmetric with respect 
to S = 7r by tuning So, e.g. Sq = t/3 for the case of 
a = 1/3. Thus the spectrum for fixed S is still of dou- 
ble degeneracy with each spin of the midgap states not 
localized at the same side of the impurity but at the op- 
posite side separatively and simultaneously. Fine-tuning 
the both phases, we may get a rich physics. Here we 
mainly demonstrate the physics of the second case with 
So = 7r/3 for the case of a = 1/3 to make the whole 
spectrum mirror-symmetric. 

Within the same harmonic trap independent of spins 
in cold atom experiments as the spinless one, we can also 
measure the LDOS for each spin using the SRRFS. As 
the relative phase S changes, the positions of spin-up and 
spin-down midgap states will be exchanged (Fig. [5]). 
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FIG. 5: LDOS for spin-up and spin-down fermions with a 
single impurity of a large attractive potential (F = —10) 
within a slowly varying harmonic trap (Vh = 0.0005). The 
midgap states for opposite spins in the bulk gaps are ex- 
changed from one side of the impurity to the other by tun- 
ing the phase shift, (a) 8 — it/3, (b) 8 — n. Here we take 
a = l/3,t=l,i = 201,V = l. 



V. PHOTONIC CRYSTALS AND THE 
ADIABATIC PUMPING 

Due to the realization of the superlattices in photonic 
crystals^, another way to realize and detect the topology 
of superlattices is with the help of photonic crystals. In 
the experiment of Ref£, the superlattice can be imple- 
mented by a set of coupled single-mode waveguides and 
the propagating light can hop between the neighboring 
waveguides. The on-site and the hopping terms of the 
superlattice can be controlled by the refraction index of 
the waveguides and the spacing between them, respec- 
tively. The propagating direction can be used to detect 
the evolution of the injected light. In this setup, we do 
not need harmonic trap and it is more easily to add a 
single on-site impurity by just tuning one waveguide's re- 
fraction index. So the above phenomena for the midgap 
states can be also seen in this photonic crystal setup. If 
we inject a light into the waveguide nearest to the im- 
purity waveguide, we will see a localized or an extended 
light signal coming out depending on the 8 we choose for 
the superlattice. To realize the spinful superlattices, two 
kinds of injected light to feel the difference of the lattice 
potential are needed. 

Besides the detection of the midgap states, we can also 
use the photonic crystals to realize the adiabatic pump- 
ing of the midgap states from one side of the impurity to 
the other. This process should rely on the 'off-diagonal' 
version of the on-site modulated superlattices, which has 
the similar properties as the on-site one and can be re- 
alized by varying the spacings of waveguides along the 
propagating direction. The Hamiltonian is as follows: 

H = + ti)clci +1 + H.c] + Fh , (3) 



where i; = Vcos(27rm + S) and t is set to the unit of 
energy. To realize the adiabatic pumping in this "off- 
diagonal" version of superlattices, the propagating axis 
is used to mimic the slowly varying of the phase shift S 
of the whole superlattice by changing the spacings of the 
waveguides along this axis. As the injected light propa- 
gates in this set of waveguides, the 5 varies as well^. Fig. 
Oshows the spectrum of the Hamiltonian ([3]) with respect 
to 5. If we inject a light into the waveguides right-nearest 
to the impurity one at S = 0.57T in the beginning, it will 
be localized at this site, and then as the light propagates, 
the S will be changed and the light will be expanded to 
other sites at 5 = 7T, and finally it will become localized 
at the left side of the impurity at S = 1.5tt. So through 
this process we can pump the midgap states from one 
side of the impurity to the other. Likewise, for spin-1/2 
fermions we can also pump the spins from one side to the 
other or exchange different spins. This kind of manipu- 
lation may be used in quantum information technique to 
control the states or spins. 
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FIG. 6: The eigenenergies versus 8 for the 'off-diagonal' ver- 
sion of superlattices ([3]) with a single on-site impurity of a 
large (F = —10) on-site potentials under the periodic mod- 
ulation, q = 1/3. The pumping process of the midgap state 
are painted red. The insets are the density profiles show the 
initial (8 — 0.57r), middle (8 — ir), and final (8 — 1.5tt) points 
of the 8 to demonstrate the transfer from one side of the 
impurity to the other. The other parameters we choose are 
t = 1, L = 201, V = 1 under PBCs. 



VI. SUMMARY 

In summary, we show that a single on-site impurity will 
induce midgap states in the ID topologically nontrivial 
superlattices. In the limit of a large impurity potential, 
these states will be localized next to the impurity site 
and behave like the edge states under OBCs. Based on 
this character of the midgap states, we propose an easier 
and more realistic scheme to detect the topology of ID 
superlattices. In cold atoms, we calculate the LDOS for 
the SRRFS measurement to demonstrate the localization 
and the transfer of the midgap states for the spinless and 
spin-1/2 ID superlattices. Alternatively, the photonic 
crystal scheme is also proposed to demonstrate the fea- 
ture of the midgap states and can be used to realize the 
adiabatic pumping. Our proposals are supplying a good 
way to show the topology of ID superlattices and may 
be also helpful to the manipulation of the midgap states 
for quantum information. 
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